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Si.  Introduction 


It  is  our  belief  that  a  thorough  understanding  of  the  mathematical  underpinnings  of  elasticity 
is  crucial  to  its  analytical  and  numerical  implementation.  For  example,  in  the  analysis  of  rotating 
structures,  if  one  attempts  to  couple  geometrically  inexact  models,  obtained  by  linearization  or  other 
approximations  to  routing  rigid  bodies,  one  can  easily  get  serious  artificial  "softening”  effects  that 
can  significantly  alter  numerical  results;  see  Simo  and  Vu-Quoc  [1986c]  for  a  discussion  (compare 
equations  (3)  arid  (5)  of  that  paper).  In  this  paper,  we  consider  geometrically  exact  models,  such  as 
the  Kirchhoff-Love-Reissner-Antman  model  for  rods  and  its  counterpart  for  plates  and  shells. 
These  models  take  into  account  shear  and  torsion  as  well  as  the  usual  bending  effects  in  traditional 
rod  and  plate  models.  Our  purpose  is  to  systematically  develop  the  Hamiltonian  structure  for  the 
dynamics  of  these  models  in  the  convective  representation.  The  convective  representation  is  chosen 
fpr  its  computational  convenience  and  for  our  planned  coupling  of  these  models  to  the  dynamics  of 
rigid  body  motion,  as  in  Krishnaprasad  and  Marsden  [19Mf 

One  of  the  topics  that  is  of  importance  in  the  foundations  of  elasticity  is  a  geometric 
formulation  of  the  equations  in  Hamiltonian  form.  JTus  form  is  extremely  useful  in  the  dynamical 
analysis  of  systems;  for  example  in  the  study  of  nonlinear  stability  (see  Holm,  Marsden,  Ratiu,  and 
Weinstein  [1985],  Krishnaprasad  [1985],  and  Lewis,  Marsden,  and  Ratiu  [1986a]),  in  bifurcation 
theory  (see  Golubitsky  and  Stewart  [1986]  and  Lewis*,  Marsden  and  Ratiu  [1986b])  and  in  the 
study  of  chaotic  solutions  (see  Holmes  and  Marsden  [1983]  and  Guckenheimer  and  Holmes 
[1983]).  Our  own  motivation  is  to  provide  additional  insight  for  work  on  routing  structures  using 
geometrically  exact  models  (see  Krishnaprasad  and  Marsden  [1986]  and  Krishnaprasad,  Marsden, 
and  Simo  [1987]).  These  motivations  are  of  course  in  addition  to  the  fart  that  these  Hamiltonian 
structures  are  of  intrinsic  interest  in  themselves  for  the  mathematical  foundations  of  elasticity 
theory.  ^ 

The  Hamiltonian  structure  for  the  material  (or  Lagrangian)  representation  of  elasticity  is  given 
in  terms  of  canonically  conjugate  variables  *  namely  the  placement  field  and  its  conjugate 
momentum,  the  momentum  density.  This  is  a  fairly  standard  result  that  is  given  in  for  example, 
Marsden  and  Hughes  [1983,  Chapter  5].  The  relation  beween  this  and  other  structures  in  spatial 
and  body  representations  is  an  important  development  that  goes  back  to  Arnold  [1966]  and  was 
developed  by  Marsden  and  Weinstein  [1974,  1982],  and  others.  A  non  canonical  Hamiltonian 
structure  for  elasticity  that  is  partially  a  spatial  represenution  is  given  in  Holm  and  Kuperschmit 
[1983],  and  a  Hamiltonian  structure  for  isotropic  elasticity  in  spatial  representation  is  given  in 
Marsden,  Ratiu,  and  Weinstein  [1984a,b].  We  concentrate  in  this  paper  on  the  convective 
representation  and  also  develop  a  Hamiltonian  formalism  for  rods  and  plates.  This  is  done  for 
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§1.  Introduction 


3 


several  reasons.  First  of  all,  no  assumption  of  isotropy  is  rr>k.ued  Second,  the  convective 
representation  is  convenient  for  the  analysis  of  coupled  nfid  body  -  beam  and  plate  structures 
Finally,  the  convective  representation  is  very  convenient  for  numerical  computations  for  | 

geometrically  exact  models  (see  Simo  and  Vu-Quoc  [1986a.b]).  A  geometric  setting  that  is  useful 
for  understanding  the  general  relation  between  the  material,  inverse  material,  spatial,  and  convective 
representations  is  given  in  Holm,  Marsden,  and  Ratiu  (1986}. 

As  we  have  mentioned,  the  noncanonical  brackets  found  in  this  paper  are  obtained  by  the  i 

general  methods  of  reduction  from  the  canonical  structure  in  material  representation,  as  in  Arnold 
[1966]  and  Marsden  and  Weinstein  [1982].  When  these  procedures  are  done  for  fixed  boundary 
problems,  one  obtains  Lie-Poisson  brackets  associated  with  the  dual  of  a  Lie  algebra  of  a 
semi-direct  product.  (See  Marsden,  Weinstein,  Radu,  Schmit,  and  Spencer  [1983]  for  a  general 
introduction  to  this  geometric  theory.)  These  sorts  of  brackets  appear  for  example  in  the  equations 
for  compressible  fluids  and  MHD  (see  Marsden  [1982],  Holm  and  Kupcr  schmit  [1983]  and 
Marsden,  Radu.  and  Weinstein  [I984a,b]).  The  geometric  reason  for  this  is  that  when  one  starts 
with  a  configuration  space  that  is  a  group  G,  and  one  reduces  the  phase  space  T*G  by  a  subgroup 
Ga  that  is  the  isotropy  subgroup  for  a  representation  of  G  on  a  vector  space  V,  then  the  resulting 
space  is  essentially  the  dual  of  the  Lie  algebra  of  the  semi-direct  product  G  ®  V.  (This  result,  due 
to  Ratiu,  Guillemin  and  Sternberg,  is  proved  in  a  sharpened  version  in  Marsden,  Ratiu,  and 
Weinstein  [1984a]  to  which  we  refer  for  the  original  references).  When  free  boundaries  are  present, 
however,  the  brackets  are  only  partially  of  the  Lie-Poisson  type.  The  geometric  setting  for  these  is 
the  "gauged  Lie-Poisson"  context  of  Montgomery,  Marsden,  and  Ratiu  [1984].  This  was  applied  to 
free  boundary  fluid  problems  in  Lewis,  Marsden,  Montgomery,  and  Ratiu  [1986].  In  this  paper, 
we  shall  not  require  the  fairly  sophisticated  context  of  the  gauged  Lie-Poisson  structures,  but  rather 
we  shall  obtain  the  results  by  a  direct  calculation.  However,  we  do  note  that  when  no  boundaries 
are  present,  the  Poisson  brackets  we  get  for  three  dimensional  elasticity  do  reduce  to  Lie-Poisson 
brackets  for  a  semi  direct  product  For  rods  and  plates,  the  brackets  also  reduce  to  the  Lie-Poisson  (  *  eopy  J 
type  in  the  cases  that  the  configuration  space  reduces  to  a  group;  for  example  this  happens  for 
torsional  motion  of  a  rod. 

The  geometric  point  of  view  adopted  in  this  paper  has  proven  particularly  useful  in  the--— - - 

numerical  solution  of  initial  boundary  value  problems.  For  the  geometrically  exact  rod  model.  for«  w 

instance,  exact  update  procedures  for  the  configuration,  stress  resultants  and  stress  couples  can  be  ^  □ 

developed  by  employing  discrete  algorithmic  counterparts  of  the  exponential  map  and  parallel  tlon _ 

transport  (see  Simo  &  Vu-Quoc  [1986a,b]).  This  results  in  algorithms  that  exactly  preserve  the - 

fundamental  physical  requirement  of  material  frame  indifference.  Similarly,  for  three  dimensional  ^^ 
nonlinear  viscoelastic  solids,  by  exploiting  the  convective  representation,  one  can  develop  lllLy  Co(SiS~ 

uncondionally  stable  and  second  order  accurate  algorithms  which  exactly  preserve  covariance  of  theu  'and/or 

Dlst  SpaelaJ. 
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continuum  formulation  (see  Simo  (19861).  Thus,  these  algorithms  go  beyond  the  notion  of 
incremental  objectivity,  as  proposed  by  Hughes  and  Winget  [1980].  Finally,  we  believe  that  the 
Hamiltonian  structures  developed  in  this  paper  will  play  a  central  role  in  future  developement, 
design,  and  stability  analysis  of  time*  stepping  integration  algorithms  for  nonlinear  elastodynamics 
which  ensure  not  only  conservation  of  energy,  (as  in  Chorin  et  aL  [19781  or  Hughes,  Liu  & 
Caughy  [1980]),  but  exactly  preserve  other  fundamental  integrals  of  motion  such  as  global  angular 
momentum. 

Acknowledgements.  We  thank  David  Fox,  Darryl  Holm,  Tom  Hughes,  Debbie  Lewis, 
Tudor  Ratiu,  and  Loc  Vu-Quoc  for  their  comments  and  joint  discussions. 


§2.  Summary  of  Three  Dimensional  Elasticity 


We  summarize  the  notation  to  be  used  in  the  description  of  three  dimensional 
elastodynamics,  following  to  a  large  extent  the  useage  of  Marsden  and  Hughes  [1983]. 


The  Configuration  Space 

Let  (2,  G)  and  {S,  Q)  be  two  smooth  Riemannian  manifolds  canying  metrics  6  and  g 
respectively.  Typically  we  have  2  C  5  and  in  fact  S  =  R^,  Euclidean  three-space  with  the 
standard  Euclidean  metric.  One  refers  to  2  as  the  reference  configuration  with  points 
denoted  by  X  €  2,  and  one  speaks  of  S  as  the  ambient  space  in  which  the  body  2  moves. 
Points  in  S  are  denoted  X  e  J.  We  shall  consider  coordinate  charts  !RA  :  2  — *  R  and  Xa  : 
S - ♦  R  so  that  the  local  coordinates  of  the  points  X  and  X  are  denoted 

XA  =  A(X)  .  for  X  €  2  and  Xa  =  Xa(x)  ,  for  X  €  S  (2.1) 

The  configuration  space  C  is  the  set  of  (orientation  preserving)  embeddings  <P  :  2 - ►  S: 

we  write 


C  =  Emb(2,  S )  (2.2) 

and  call  the  set  <P(2)  the  current  confuguration.  It  is  known  that,  when  suitably  tcpologized, 
C  is  a  smooth  infinite  dimensional  manifold  (see  Abraham,  Marsden,  and  Ratiu  [1983],  Ebin  and 
Manden  [1970]  and  references  therein). 

The  tangent  space  to  C  at  a  configuration  <P  is  given  by 

TfC  =  {  V  :  T2 — ►  T5  |  V(X)  e  T^(x)5forall  X  €  2)  (2.3) 


Kinematics 

A  motion  is  a  curve  of  configurations;  we  let  <Pt  be  the  configuration  at  time  t  and  write 
<pt(X)  -  <P(X,t).  Given  a  motion  <Pt,  we  define  the  following  quantities: 

(i)  material  velocity  :  Vt  c  T^C  given  by 
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Vt(X)  :=  (PCX.  t)  (2.4a) 

(ii)  spatial  velocity  :  Vt  €  3C(<Pt(®))  [the  space  of  vector  fields  on  tpt((B)]  is 
defined  by 

vt  =  Vt  .  V'  (2.4b) 

(ill)  convective  velocity  :  e  X  [fB)  is  defined  by 

Vt  =  <P;(Vt)  ;=  T<Pt_1  •  Vt  •  <Pt  =  T^"1  •  Vt  (2.4c) 

The  deformation  gradient,  denoted  Ft  is  defined  to  be  the  tangent  map  of  (Pt;  we  write 
Ft  =  T<Pt.  In  coordinates. 


F* 


Proposition  2.1.  The  convected  velocity  is  minus  the  spatial  velocity  of  the  inverse 
motion  (P^1  :  S - *  ®  ;  I.C., 


^t 


<Pt 


Proof.  Applying  the  chain  rule  to  the  identity  X  s  <P_1  (<P(X.  t),t)  gives 

if,'  -1 

ST' 9,  ‘Tf,  ’•  V,  =  0 

and  so  the  result  follows  by  noting  that  =  T(Pt”*  •  Vt  •  <Pt.  ■ 


(2.5) 


The  Metric  and  Convected  Metric  Tensors 


The  convected  metric  tensor  is  defined  by  the  pull-back  relation: 


W  L^l**  W  5W  lNT.Hfr 


W  IAJ  \mj  y^,\  { ^ 
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One  refers  to  as  the  right  Cauchy-Green  tensor.  In  coordinates, 


CAB  =  F*AFbB9ab*^- 


(2.6b) 


Associated  with  the  convected  metric  Ct  we  define  a  symmetric  Levi-Cevila  connection  V 
to  be  the  connection  of  the  metric  Ct;  i.e.  by  the  standard  relation  for  the  Christoffel  symbols: 


-  VA-W*  A  ♦  rADBV A'H'0 

where 

r ,o0  =  cBEr ^DE 

and 


1  F  dCAC  ^CBC  ^CAB  I 

2  L  exB  dxA  dxc  J 


For  the  case  of  Euclidean  space  one  has  the  formula 


d2<P8 

dXAdX° 


(F-1)B 


a 


(2.7a) 


(2.7b) 


(2.7c) 


(2.8) 


(See  Marsden  and  Hughes  [1983],  p.  31.) 


Acceleration  Vectors 

In  addition,  associated  with  the  motion  <Pt  ,  we  define  the  material  acceleration  At 
and  the  spatial  acceleration  Bt  by  the  expressions 

At  =  d2<pt/et2  =  evt/dt  and  8t  =  At  •  tp'1  (2.9) 

The  convected  acceleration  is  defined  by  the  pull-back  relation 


(2.10) 


k*i*«*£j 


*  *  *->  ».  *  v  A  ;■»  /V  V  *  VT»  V“J»^ 
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Proposition  2.2.  The  connected  velocity  and  acceleration  are  related  by  the  formula 


*vy, 


(2.11a) 


In  coordinates,  A,  is  given  by 

j^sd^/dt  *  ♦  J"cDA'k'cV°  (2.11b) 

Proof.  Recall  (Marsden  and  Hughes  [1983]  p.  33)  that  the  spatial  velocity  and  accelerations  are 
related  by 

dv, 

»t  =  -5T  *  V*,v>  •  (212) 


Now  pull  back  the  relation  (2.12)  by  <Pt;  one  gets 

=  vtv**,v«  =  vvvf 

and  by  the  Lie  derivative  formula  (Marsden  and  Hughes  [1983],  §1-6), 


Adding  (2.12)  and  (2.13)  gives  (2.11).  ■ 

The  Stress  Tensor  and  Covariance 

We  assume  the  existence  of  a  stored  energy  function  W  :  *  C*  My  — »  R.  where 

is  the  space  of  Riemannian  metrics  on  S,  and  is  the  space  of  Rietnannian  metrics  on 
%  of  the  form 


W  =  W(g,  F.G) 


(2.14) 


where  W  depends  only  on  the  point  values  of  Q,  F,  and  G.  This  is  in  keeping  with  the 
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classical  assumption  that  the  stored  energy  function  in  an  elastic  material  depends  on  the 
configuration  <P  only  locally  through  the  point  values  of  the  deformation  gradient  F  (see  Marsden 
and  Hughes  (1983],  $$3.2  and  3  .3).  The  dependence  of  the  stored  energy  W  on  the  metric  tensor 
Q  is  essential  to  introduce  the  notion  of  covariance.  Consider  an  arbitrary  superposed  spatial 
diffeomorphism  T\  :  S  — *  S',  then  we  say  that  W  is  covariant  if 

W(g,  F,G)*n  =  W(T\«g.  Tt(*F,G)  (2.14) 

This  assumption  implies  that  W  depends  only  on  the  point  values  of  the  Cauchy-Green  tensor  C. 
That  is,  there  is  a  function  of  the  point  values  of  metrics  on  which  we  shall  denote  by  W  such 
that 


W(g.  F.G)  =  W(C,G)  (2.15) 

Let  CJ  be  the  Cauchy  stress  tensor  and  let  Z  ~  Jac (<P)0.  (where  Jac(<P)  is  the 
Jacobian  determinant  of  <P),  be  the  Kirchhoff  stress  tensor.  Finally,  the  symmetric 
Piola-Kirchhoff  stress  tensor  is  defined  by  the  pull-back  relation 

S  =  <p*(r)  =  Jac(tp)«p**(cj)  (2.16) 


We  then  have  the  constitutive  relations 


z  =  2p 


aw(g,  F,G) 

nf  dg 


and 


„  aWc.s) 

S  =  2P«,5c - 


(2.17) 


where  pref  is  the  density  in  the  reference  configuration.  The  relation  (2. 1 7 1 )  is  refeiTed  to  as  the 
spatial  Doyle-Ericksen  formula  (see  Marsden  and  Hughes  (1983]  $3.3,  and,  for  the  material 
counterpart,  Simo  and  Marsden  (1984]).  In  terms  of  the  Lagrangian  strain  tensor  E  =  (C 
-  G)/2.  and  the  Euterian  strain  tensor  6  ~  <P«(E)  =  (Q  -  b"1)/2,  where  b-1  = 
<P»(G)  is  the  Finger  deformation  tensor,  formulae  (2.17)  read 


Z  =  p 


dW(e.F.G) 


ref 


de 


and 


S  =  Pref 


dW(E.G) 

eE 


(2.18) 


Note  that  the  dependence  of  W  on  the  material  metric  tenso  G  has  been  explicitly  assumed  in 
equation  (2.14)  but  that  G  is  treated  as  a  parameter  as  far  as  the  covariance  assumption  is 
concerned. 
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Next,  we  consider  the  invariance  group  of  W  on  the  left  and  on  the  right.  We  observe  that 
the  form  of  the  stored  energy  function  W  =  W(fl,  F,  6)  is  left  invariant  under  the  action  of  the 
group  Diff(j)  of  all  diffeomotphisms  of  S  onto  5.  Indeed,  this  follows  directly  from  the 
covariance  assumption.  Thus,  we  regard  the  function  W(C,  6)  as  being  obtained  by  left 
reduction  from  the  function  W.  This  relation  is  key  to  the  derivation  of  the  Hamiltonian  structure 
in  the  convective  representation  and  will  be  explored  in  detail  in  the  next  section. 

The  Convective  Hamiltonian 

We  conclude  this  section  with  a  discussion  of  the  Hamiltonian  in  the  convected  picture.  In 
the  Hamiltonian  formalism  for  the  material  picture,  the  kinetic  and  potential  energy  are 
expressed  in  the  variables  on  T *C:  i.e.,  in  the  variables  $  e  C  and  its  conjugate  momentum 
€  T*C  which  is  related  to  the  material  velocity  by  the  coordinate  formula 

Ma(X)  =  gab(x)Vb(X)P|t£(X)  d3X  (2.19) 

so  that  M,j>  is  a  one  form  density  covering  <P.  We  obtain  the  expressions  for  the  kinetic 
energy  in  the  convective  representation  by  simply  changing  variables  as  follows.  The  square  length 
of  the  material  velocity  is  given  by 

|Vt(X)|*  =  V8(X)Vb(X)gab(x)  (where  X  =  <J>t(X)) 

--  F8AFbB1/AvBgab(x) 

e  Cab0/A'V/B.  (2-2°) 

Introducing  the  inner  product 

<P,W>  ::  f  P„fCABVAVSd3X  (2.21) 

3 

on  the  space  of  convective  velocity  fields,  the  kinetic  energy  may  be  written  as 

-  \<v.v> 


(2.22) 
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This  kinetic  energy  induces  a  function  on  the  cotangent  bundle  in  a  standard  manner.  The 
corresponding  momentum  induced  by  this  kinetic  energy  via  the  Legendre  transform  is  given  by 

n  :=  Piti^^X  (2.23) 

where  ^  denotes  the  index  lowering  operation  so  that  yb 

is  given  in  coordinates  by  7^  = 

and  d3X  is  the  coordinate  volume  element  on  H.  Similarly,  following  standard  notation 
(see,  for  example,  Abraham,  Marsden,  and  Ratiu  [1983]),  the  index  raising  operation  is  denoted  by 
*,  so  that  is  given  in  coordinates  by  =  CABlMg,  where  CAB  denotes  the  inverse  matrix 
to  CAB.  Substitution  of  (2.23)  into  (2.22)  gives 

x  =  J  f  2  J~  !tfA*4CABd3X  (2.24) 

Note  that  K  is  now  a  function  of  Tl  and  C  alone.  We  also  note  that  the  Irinetic  energy  is  just  one 
half  the  square  length  of  the  momentum  density  in  the  metric  on  the  space  of  convective  momentum 
densities  that  is  induced  by  (2.21).  Thus  the  total  Hamiltonian  is  given  by  the  expression 

H(n.  C.  6)  =  M  7T  ^A*4CABd3X  ♦  f  Preftf(C.  6)  d3X  (2.25) 
2  Pref  4  'B 

We  think  of  this  form  of  the  energy  as  being  the  energy  induced  on  the  original  space  of  M, 
g's  after  quotienting  by  the  group  of  spatial  diffeomorphisms  Diff(5)  (and  again  G  enters 
parametrically).  This  idea  is  central  to  the  reduction  procedure  that  will  be  explained  in  the  next 
section. 

Convected  Equations  of  Motion 

The  convected  equations  are  obtained  by  either  pulling  back  the  equations  of  motion  in  spatial 
coordinates  or  by  a  coordinate  calculation  of  dfl/dt ,  etc.  From  balance  of  linear  momentum,  one 
gets  the  following: 

n  =  C-DivcS  ♦  [WYYl  (2.26a) 

whereas  balance  of  angular  momentum  yields  the  standard  symmetry  condition  S  =  S^.  Finally 
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one  gets  the  additional  relation: 

t  -  LyC  (126b) 

Here,  Ly  C  denotes  the  Lie  derivative  along  the  flow  of  the  converted  velocity  field 
Equation  (2.26b)  follows  directly  from  the  relation  between  flows  and  Lie  derivatives  and  the 
relation  Ct  =  Another  argument  is  as  follows.  Recalling  (2.4c)  and  (2.6a),  we  have 

LpC-Lf.(v)9"(C)  =  <P"(  Lvg)  (2.27) 

Since  Lvg  =  2d,  where  d  is  the  rate  of  deformation  tensor,  (2.26b)  follows  from  (2.27)  by 
recalling  the  well  known  relation 

c  =  2<p;(d). 

It  is  equations  (2.26a)  and  (2.26b)  that  we  shall  show  are  Hamiltonian  in  the  next  section. 


« J « U  MSJ  WM 


aiijuu iijA,l^|,j„.n  mj,  i, 


$3.  Hamiltonian  Structure  of  Three  Dimensional 
Elasticity  in  the  Convective  Representation 

In  this  section  we  show  that  the  equations  of  elastodynamics  in  the  convective  representation 
are  Hamiltonian  relative  to  a  non-canonical  Poisson  structure  on  the  space  of  pairs  [TV,  C), 
where  n  is  the  convective  momentum  density  and  C  is  the  Cauchy-Green  tensor,  as  in  the 
preceeding  section.  This  means  that  the  equations  of  elastodynamics  are  equivalent  to  this 
condition:  for  any  function  f(TT,  C), 

f  =  {f.  H }.  (3.1) 

where  H  is  the  Hamiltonian,  given  by  equation  (2.25).  The  bracket  (  )  satisfies  the  usual 
conditions  for  a  Poisson  bracket,  including  Jacobi's  identity  (see,  for  example,  Marsden  et  al. 
[1983]  for  some  generalities). 

The  bracket  (,}  appearing  in  (3.1)  will  be  obtained  by  reducing  the  canonical 
bracket  on  T *C  by  the  group  of  spatial  isometries  of  the  metric  fl  on  S  .  Equivalently,  as 
in  Marsden,  Ratiu,  and  Weinstein  [1984a,b],  we  can  add  the  metric  Q  as  a  parameter  and  reduce 
T*(C  *  iWj)  by  the  left  action  of  Dlff(i).  This  reduction  procedure  will  be  explicitly 
explained  by  direct  calculation  in  what  follows.  Before  reading  this  section,  the  reader  may  find  it 
helpful  to  first  review  the  parallel  case  of  rigid  body  dynamics  in  $4. 


The  Canonical  Bracket  on  the  Material  Phase  Space 

We  start  with  the  canonical  bracket  on  T*C.  the  space  of  deformations  <P  and  their 
canonically  conjugate  momenta  the  material  momentum  densities.  Now  let  f  and  g  be 
functions  of  the  pairs  (M,  C)  .  Define  the  new  functions  f  .  g  on  1*C  by 

7 (<f>.  M*,)  =  <P*fl)  (3.2) 


where 


<P*M,  =:  n  =  (T«P)t  • 


(3.3) 


Here.  (T<?)t  =  is  the  transpose  of  the  deformation  gradient  and 
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<p*g  =  c  =  ftf 

as  before.  In  coordinates,  (3.3)  reads 


(3.4) 


itfA(X)  =  F8A(X)Ma(X)  (3.3) 

and 

CAB(X)  =  F'a()OF6b(X)9,i,(x).  (3.4) 


where  X  =  <P(X)  .  We  define  g  in  a  similar  way.  Now  we  form  the  canonical  bracket  on 

TV 


(TTl  -  f  (iI--iS-.il- JLLW 

1  ,8'  "Ja'sip  sm„  s«  sm,,'0* 


where  Sf  /5<P  is  the  functional  derivative.  (This  Poisson  bracket  is  written  as  if  <P  and 
were  independent  variables;  this  requires  a  little  caution  since  T*C  is  not  simply  a  product  space. 
However,  it  is  well  known  how  to  deal  with  this  situation;  see  Abraham  and  Marsden  [1978]  and 
Marsden  and  Hughes  [1983]).  The  bracket  (3.5)  is  now  to  be  expanded  using  the  chain  rule  and 
the  result  expressed  in  terms  of  f.Q.TV,  and  C.  This  will  produce  the  desired  reduced 
bracket.  General  reduction  theory  (Marsden  et  al.  [19831,  Marsden  and  Ratiu  [1986]  and 
references  therein)  shows  that  the  elastodynamic  equations,  which  are  Hamiltonian  on  T  C  (or 
Lagrangian  on  TC),  reduce  to  Hamiltonian  equations  (3.1)  on  our  convective  phase  space  T, 
the  space  of  (TT,  C)  's. 

Before  proceeding  with  this  calculation,  we  make  some  additional  remarks.  First,  if  C  were 
a  group,  we  would  expect  the  reduced  bracket  to  be  of  a  special  type,  namely  a  Lie-Poisson  bracket 
on  the  dual  of  the  Lie  algebra  of  a  semi-direct  product;  see  Marsden,  Ratiu,  and  Weinstein 
[1984a,b].  We  shall  see  analogues  of  such  a  structure  here.  Because  the  boundary  of  can 
move,  C  is  not  a  group  and  indeed  our  bracket  differs  from  a  semi-direct  product  Lie-Poisson 
bracket  only  in  boundary  integral  terms.  For  another  bracket  in  spatial  representation,  see  Lewis, 
Marsden,  Montgomery,  and  Ratiu  [1986];  see  also  Holm,  Marsden,  and  Ratiu  [1986].)  The 
boundary  conditions  here  can  either  be  displacement,  traction,  or  a  combination.  In  the  former 
case  we  restrict  C  ,  which  imposes  corresponding  restrictions  on  V  and  TV;  specifically, 
if  a  portion  of  dS  is  fixed,  TV  will  vanish  on  that  portion.  For  traction  boundary  conditions 
we  add  a  corresponding  boundary  term  to  the  Hamiltonian  (as  in  Marsden  and  Hughes  [1983],  for 
example).  These  boundary  terms  are  to  be  taken  into  account  in  the  calculation  of  functional 
derivatives,  as  in  Lewis,  Marsden,  Montgomery,  and  Ratiu  [1986]). 

A  second  remark  is  this.  The  bracket  in  purely  spatial  terms  seems  to  be  possible  only  for 
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isotropic  elasticity;  see  Marsden,  Ratiu,  and  Weinstein  [1984a,b].  No  such  restriction  is  needed 
for  the  convective  representation.  The  reason  for  this  discrepancy  is  simply  that  material  frame 
indifference  gives  us  a  covariance  symmetry  with  which  to  reduce  by  Diff(5)  .  The 
corresponding  reduction  by  D  i  f  f  {$)  requires  isotropy  of  the  material. 

A  third  remark  is  that  given  a  dynamic  solution  M(X,  t),  C(X,  t),  the  original  motion 
<P(X ,  t)  can  be  reconstructed  by  constructing  V  using  (2.23)  and  then  solving  the  ordinary 
differential  equation  for  Proposition  2.1: 

S  (3.6a) 

and  letting 

<Pt  =  (3.6b) 

This  is  a  special  case  of  the  general  reconstruction  procedure  that  is  used  in  reduction  theory.  (See 
Abraham  and  Marsden  [1978]). 

The  reduced  Poisson  structure  is  given  in  the  following  theorem.  The  notation  is  as  follows: 

(i )  S  f/STt  is  the  vector  field  on  OB  such  that  for  all  variations  571  (which  are  also  a 
one  form  density,  like  71), 

Df-Sn  :=  f(71  ♦  tSTt,  C)  1  =  f  4Sf-SHd3X.  (3.7a) 

<36  t  e=0  J  <3  brt 

(ii)  6 f /SC  s  the  symmetric  two  tensor  (indices  up)  density  on  CB  such  that 

Df-SC  :=  4-  un.  C  ♦  ESC)  I  =  f  4L  !  6C  d3X  (3.7b) 
oe  I  e=o  J 

for  all  variations  8C  (which  is  a  symmetric  two  tensor,  with  indices  down,  on  ®.) 

(iii)  0 1  V(*)  denotes  the  divergence  relative  to  the  coordinate  chan  {XA}  with  associated 
(reference  configuration)  metric  G ^(X),  whereas  Divc(0  denotes  the  divergence  operator 
relative  to  the  convected  metric  C^g(X).  In  addition,  covariant  differentiation  relative  to  the 
convected  metric  C^g  is  denoted  in  the  remainder  of  this  section  simply  by  V  or  in  coordinates 
with  a  vertical  bar,  for  instance,  DivcV  =  T*  |  A  =  tr[Vl>]. 

Let  us  denote  by  10*  the  skew-symmetric  part  of  a  second  rank  tensor  10  .  Then  w  e 
define  the  convected  vorticity  W  by 
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W  *  H  s  i^n^H  ,  for  all  vectors  H  e  R5  . 

2 


(3.8) 


Theorem  3.1  The  reduced  Poisson  bracket  on  IP  is  given  by 


where  N  is  the  outward  unit  normal  on  d®,  L  is  the  Lie  derivative  (recall  that  Sg/SC  is  a  two 
tensor  density  with  indices  up),  and  where  ( ,  ]  denotes  the  usual  Lie  bracket  of  vector  fields. 

Note  that  for  displacement  boundary  conditions.  (Sf/6H)'N  =  0  ,  and  so  the  two 
boundary  terms  in  (3.9b)  disappear,  resulting  in  a  Lie-  Poisson  bracket  for  the  semi-direct  product 


-  *-•  *-  *  y',yj 


1  ”v  ’  jv.  nr 
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X(®)  ©  (2-tensor densities). 

We  turn  now  to  the  proof  of  the  theorem.  For  simplicity  of  notation,  we  just  write  f  for  f . 
Lemma  3.2.  The  following  formulae  hold: 


2i_  -  .Q  rbJ22L.  a/A  21-1 

5(pa  9ab  A*-  dCAB  J 


2L  =  FaA  2L 

8Ma  a  e*fA  * 


(3.10a) 


(3.10b) 


Proof.  The  second  formula  follows  at  once  from  the  chain  rule.  To  prove  the  first  one  we 
proceed  as  follows.  Denoting  by  <  •  ,  •  >  the  standard  l^-pairing,  we  have 


d*fA  .A 


5  t  dv 


(3.11) 


wh^re  use  has  been  made  of  the  relation 


M.  =  (F"’)c.SWb  =  (F''  )C.CCA«*  =  9.bFbA^  •  (3  12) 


Next,  recall  that  if  P  and  S  are  related  through  P,A  =  F#bSBa  ,  then  one  has 


D1V  P  =  T^*DivcS  . 

Using  the  divergence  theorem  (and  assuming  S*?*  is  zero  on  d®  )  we  have 

<T=-W,>  5 


(3.13) 


(3.14) 


\  •. 
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Hence 


er 


14-  =  -ab-F*Divc[2 

6<p  *  CL  ecb 


♦  n*®  —  ] 

dtt  J 


which  completes  the  proof.  ■ 

Making  use  of  this  proposition,  and  the  identity  |  g  ■  0,  we  have 


/ il  iiN  f 

’  6Ma'  ^  *B 


Thus 


Since  6  :  A  s  6  :  [A3A  for  any  A  e  GL(R3)  (see  equation  (3.8)),  and  § 
(the  skew  symmetric  3*3  matrices),  we  have 


vn 


9 


df  df 

“5n  5FT 


9 


U> 


r  df  dg  1 

l-Zn'Sn* 


(3.15) 

(3.16) 


(3.17) 

e  so(3), 

A 

iil 

dnJ 

(3.18) 


where  A  :  R3  — *  so(3)  is  the  usual  isomorphism  (recalled  explicitly  in  section  4). 
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Substituting  (3.18)  into  (3.17)  gives  (3.9a).  Formula  (3.9b)  is  obtained  from  (3.9a)  by 
integration  by  pans  (the  divergence  theorem).  ■ 


Corollary  3  J .  Hamilton’s  equations  f  =  { f .  H }  are  equivalent  to  the  connected  equations 


Pref< &  *VVV)  =  DivcS 

(3.19a) 

[or 

n  =  C-Divcs  ♦  ivv]Jn 

(3.19a)] 

and 

t  =  LyC 

(3.19b) 

This  follows,  as  was  remarked  earlier,  by  the  general  theory  of  reduction.  It  also  can  be 
checked  by  a  direct  calculation. 
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In  what  follows  we  summarize  some  basic  notation  and  elementary  properties  of  the  rotation 
group  needed  for  subsequent  developments.  For  a  more  detailed  account  we  refer  the  reader  to 
standard  textbooks,  such  as  Abraham  and  Marsden  [1982,  $4.1]  and  Choquet  and  DeWitt-Morette 
[1982,  pp.  181*194].  We  then  review  the  Hamiltonian  structure  of  rigid  body  dynamics. 


The  Rotation  Group  and  its  Lie  Algebra 


Following  standard  usage,  we  denote  by  S0(3)  the  Lie  group  of  proper  orthogonal 
transformations,  i.e., 


SO(3)  :=  {A  :  R3  — »R3  A  is  linear,  ATA  =  I  ,  and  det  A  =  1  };  (4.1) 


S0(3)  is  a  compact  subgroup  of  the  general  linear  group  GL(3)  .  Its  Lie  algebra  is  SO(3) 
:s  (£  :  R3  — .R3  |6  is  linear  and  6  ♦  $T  =  0),  the  set  of  all  skew- symmetric  tensors. 
The  notation  is  as  follows:  In  coordinates,  relative  to  an  orthonormal  basis  {e,}  in  R3.  we 
write  6  =  ^'j®i  ®  ©  :  ©'6,,  and,  in  matrix  notation. 


0 

-e3 

e2 

’©r 

e3 

0 

-e1 

.  {©‘}  = 

©2 

-02 

• 

e1 

0 

©3 

(4.2) 


Recall  that  So(3)  and  R3  are  related  through  the  isomorphism  :  R3  — *  so(3) 

,  defined  by  the  relation  8h  r  0  ■  h  .  for  any  h  c  R3  .  Here,  0  c  R3  is  the  axial 
vector  of  6  c  So(3)  ,  and  ■  denotes  the  ordinary  vector  product.  Physically,  A  c 
S0(3)  represents  a  finite  rotation  while  infinitesimal  rotations  are  linearized  rotations  about  the 
identity.  Geometrically,  one  speaks  of  SO(3)  as  the  tangent  space  of  S0(3)  at  the  identity  1 
c  S0(3)  ,  and  employs  the  notation  So(3)  =  T,S0(3)  . 


Left  and  Right  Representations  of  the  Tangent  Space  to  SO(3) 


Given  any  A  «  S0(3)  ,  elements  of  the  tangent  space  T^S0(3)  at  a  point 
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A  €  S0(3)  arc  represented  in  two  alternative  forms: 

(i)  Left  invariant  vector  fields  defined  by  the  relation  $A  '*  t,la6  ,  where  6  e 
So(3)  ,  and  LA  :  SO(3)  — »  SO(3)  denotes  left  translation;  i.e.,  LaAi  :=  AA^  . 
A  simple  calculation  yields  @A  =  A^  .  Accordingly,  we  have  the  following  identification  of 
TAS0(3): 


TaS0(3)  =  {@A  ;=  A§  |  for  any  §  e  So(3)}  .  (4.3a) 

Geometrically,  an  element  §A  e  TaS0(3)  corresponds  to  a  finite  rotation 
superimposed  onto  an  infinitesimal  rotation  S  e  So(3)  . 

(ii)  Right  invariant  vector  fields.  The  characterization  is  identical  to  that  in  (i),  with  left 
translations  replaced  by  right  translations  RA  :  S0(3)  — — ♦  S0(3)  defined  as  RaA^  = 
AjA  .  This  leads  to  the  representation 

TaS0(3)  :=  (0A  :=  ©A  |  for  any  §  e  So(3)}  .  (4.3b) 

Geometrically,  an  element  $Ae  TaS0(3)  represents  an  infinitesimal  rotation 
superimposed  onto  a  finite  rotation  represented  by  A  .  Following  the  general  conventions  we 
use  in  elasticity,  left  invariant  vector  fields  at  a  point  A  are  denoted  by  upper  case  letters;  i.e., 
6 A  •  We  think  of  S  €  S(>(3)  ,  with  axial  vector  9  €  R3  ,  as  a  material  object  with 
coordinates  ©*  relative  to  a  material  basis  {Ej}  .  In  the  context  of  rigid  body  dynamics  one 
often  speaks  of  body  representations  and  body  coordinates  relative  to  the  body  frame  It,} 
defined  by  t,  :=  AE,  .  On  the  ot  er  hand,  right  invariant  vector  fields  at  a  point  A  are 
denoted  by  lower  case  letters;  i.e.,  0A  .  One  thinks  of  0  e  SO(3),  with  axial  vector  0  e 
R3.  as  a  spatial  object  and  coordinates  0*  relative  to  a  spatially  fixed  basis  {e,}  .  In 
the  context  of  rigid  body  mechanics  one  speaks  of  the  spatial  representation  . 


The  Dual  of  the  Lie  Algebra 

We  also  recall  that  so(3)*  ,  the  cotangent  space  at  the  identity  1  €  S0(3)  .  is 
isomorphic  to  IR3,  *]  and  to  SO(3)  via  the  dot  product;  that  is,  we  also  identify  elements  of 
So(3)*  with  skew  symmetric  matrices  and  use  the  map  *  and  the  pairing  given  by 


I,  |U|.  I  J»| 


kk<  -  *_*  u  rj  <A  .UUn 
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ft(W)  =  j  rt:W  ■  n-w  . 


where  ft:W  =  trace  [TTtW],  and  we  have  used  the  notation  A:B  =  A^j  Bj^  for  any  A, B 
€  GL(3).  Thus,  the  duality  pairing 


<*.*>:  T*S0(3)  x  TaS0(3)  — *  R 


(4.5a) 


is  defined  by 


<HA.  *A>  ==  trace  [nAWA]  . 


(4.5b) 


Note  that  this  duality  pairing  is  left  invariant  in  the  sense  that 


<ftA.  WA>  :=  j  trace  [ftTATAft]  ■  ^  ft  :  W  =  <ft,W> 


(4.5c) 


The  Exponential  Map 


Finally,  we  recall  that  the  straight  line  E  *-»  E0  €  TjS0(3)  ,  for  E  >  0,  is  mapped  by 
the  exponential  map  onto  the  curve 


t  *-»  explte]  =  [  Z  e*  ]  e  so(3) . 


Note  that  for  the  case  of  6  e  SO(3)  one  has  the  following  explicit  fonnula  often  credited 
to  Rodrigues  (Goldstein  [1980,  p.  165]): 


exptll5  ,.^Le.  » 


This  fonnula  is  of  fundamental  importance  in  the  numerical  solution  of  initial  boundary  value 
problems  for  finite  strain  rods  (Simo  and  Vu-Quoc  [1986a]). 

Let  us  now  recall  the  Hamiltonian  structure  of  rigid  body  dynamics  in  body  (*  convective) 


I 


i 


KV 
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representation.  This  is  done  purely  for  the  reader's  convenience  to  see  the  parallel  with  the 
developments  in  the  preceding  and  subsequent  sections.  It  is  also  useful  to  keep  this  parallel  in 
mind  for  the  problem  of  coupled  dynamics  (Krishnaprasad,  Marsden,  and  Simo  [1986]).  For 
further  details  on  rigid  body  dynamics  in  this  context,  see  Marsden,  Ratiu,  and  Weinstein  [1984b], 
The  configuration  space  is  C  *  S0(3)  ,  and  the  Hamiltonian  H  :  T*S0(3)  — *  R 
is  given  by 

H  =  I  n-  l^n  .  (4.8) 

Here,  1  is  the  inertia  dyadic  defined  in  terms  of  the  Lagrangian  coordinates  X  =  (X 1 ,  X2, 

X3)  as 


I  :=  f  P„,(X)  [  |xpl  -  X  «  X  ]  d3X  ,  (4,9) 

Q 

and  TT  is  the  body  momentum  defined  as 

TT  :s  l  [AT  A  ]  .  (4.10) 

We  observe  that  H  is  invariant  under  spatial  isometries ;  that  is,  invariant  under  the  left  action  of 

$0(3)  .  Reduction  by  this  symmetry  amounts  to  considering  functions  f  :  T*S0(3) - ♦  R 

of  the  form 

f(A,  rtA)  -  7(ATrtA)  ■  7(ft)  .  (4.11) 

The  canonical  bracket  on  T"SO(3)  is  given  in  terms  of  the  duality  pairing  (  •  ,  O  : 
T^S0(3)  *  TaS0(3)  — ►  R  defined  by  (4.5)  as 


The  Lie  Poisson  bracket  on  the  reduced  space  T  ~  T"SO(3)/SO(3)  is  obtained  with  the  aid 
of  the  results  given  below 


Proposition  4.1.  The  following  formulae  hold 
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—  A 


I • 


(4.13) 


d  f 


•  *[&] 


dlt 


(4.14) 


where  f  (ft)  =  ((FT). 


Proof.  To  prove  (4.13),  observe  that  for  any  8  A  e  Ty^S0(3)  we  have  the  left  representation 
8A  =  AS6.  where  s$  e  SO (3).  Thus,  by  the  chain  rule  and  the  left  invariance  of  the 
duality  pairing  we  have 


df 


dn 


df 


=  <  -  It  .  Sfi  > 

dn 


d  f 


Since  — —  €  So(3),  it  follows  that 

dft 


<  ~ST~  -  >  =  <■  1  [n.  |£-].  S$> 


dA -  '  2 

So  that,  on  setting  f  (ft)  =  f  (TT)  and  recalling  the  Lie  bracket  relation 


AS  -  6A  =  [A  *  B3“, 


we  obtain 


<T7T  •  5A  >  =  T  A  t  n  “  “5H  1  ’  8A>- 


and  result  (4.13)  follows.  An  analogous  calculation  holds  for  the  formula  (4.14). 
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Substitution  of  (4.13)  and  (4.14)  into  (4.12),  and  use  of  standard  vector  product  identities 
yields  the  result 

(f.g)  =  -n.[AL  .  |i]  .  <4.i6> 

which  is  the  standard  Lie-Poisson  bracket  for  rigid  body  dynamics.  The  equation  f  =  {f,  H}  is 
then  easily  checked  to  be  equivalent  to  Euler's  equation  for  rigid  body  dynamics:  IT  =  Q  *  TT, 
where  fi  =  atA  is  the  body  angular  velocity. 


$5.  Geometrically  Exact  Finite  Strain  Rod 


The  static  version  of  the  rod  model  summarized  below  goes  back  essentially  to  Reissner 
[1973]  who  modified  the  classical  Kirchhoff  -  Love  model  (see  Love  [1944])  to  account  for  shear 
deformation.  An  equivalent  model,  formulated  as  a  constrained  director  theory  -  the  so-called 
special  theory  of  Cosserat  rods  -  is  due  to  Amman  [1974]  -  see  also  Antman  and  Jordan  [1975], 
Amman  and  Kenny  [1981],  and  Amman  [1984]  for  some  applications.  The  dynamic  version  along 
with  the  parametrization  discussed  below  is  given  in  Simo  [1985].  For  completeness,  a  brief 
account  is  outlined  next. 


The  Configuration  Space 

From  a  physical  standpoint,  the  configurations  of  a  rod  deforming  in  the  ambient  space  R3 
may  be  defined  by  specifying  :  (i)  The  position  of  its  line  of  centroids  by  means  of  the  map  <P  : 
[0, 1]  — ♦  R3,  and  (ii)  The  orientation  of  cross  sections  at  points  S  e  [0  .  L].  This  can  be 
done  using  the  orientation  of  a  moving  basis  { t ,(S)  |  I  s  1 ,2.3 }  attatched  to  the  cross  section 
relative  to  a  fixed  frame  { Et  |  1  =  1 ,2,3),  referred  to  as  a  material  frame  in  what  follows. 
The  moving  basis  is  described  by  means  of  an  orthogonal  transformation  A  :  [0  ,  L]  — ♦  S0(3) 
such  that  tj(S)  =  A(S)Et. 

If  we  view  the  rod  as  having  a  finite  cross  section  given  by  a  compact  set  Q  C  R2 
with  a  smooth  boundary  6Q,  then  the  placement  of  the  rod  in  its  reference  configuration  is 
determined  by  a  map  4>o(S)  =  (<Po(S),  Ao(S))  in  such  a  way  that  the  corresponding  set 
ocupied  by  the  rod  is  given  by 

2 

2  =  {X  €  R3|  X  =  <P0(S)  ♦  Z  £.aA0(S)Eo,  where  (£.'.£.2.S)  €  Q  *  [0.  L]} 

a*l 

(5.1a) 

Without  loss  of  generality  we  can  assume  that  |<Po '(S)|  =  1 ,  so  that  L  is  the  length  of  the 
reference  line  of  centroids.  Typically,  one  chooses  A(j(S)  so  that  T  j(S)  :=  Ao(S)Ej  is  the 
Frenet  frame  associated  with  the  curve  <Pq;  that  is, 

Tj  :=  <Po\  T,  :=  fo'y|9oi  and  T2  :=  T3  *  T,  (5.1b) 


According  to  these  ideas,  we  take  for  our  configuration  space  the  set 


•  ww  \r* 


- - Tt  I..OI.., 
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C~-  {♦  =  (<P.  A)  |  CO.  L]  — ►  R3  *  S0(3)  }  (5.2) 

For  simplicity,  we  shall  assume  that  and  A  satisfy  pure  displacement  boundary  conditions 
in  what  follows;  i.e.,  <P  J  5  _  q_|_  and  A  j  5  =  o,t  811  Prescribed.  Hence,  the  tangent  space  at 
the  identity  configuration  is  given  by 

T  C 

identity4' 

=  {(S<p.  s§):  [0.  U  — ■*  R3  *  so(3)  |  5<P js  .  0§L  =  0 

andS@|g  s  q  l  •  0).  (5.3) 


Left  and  right  invariant  tangent  vector  fields  at  a  configuration  ♦  e  Care  defined  below  in  the 
standard  fashion  by  employing  left  and  right  translations.  Finally,  associated  to  any  configuration 
♦  €  C  one  defines  its  arc-length  by  the  mapping 

S  e  [0  ,  L]  *-*  S  =  s(S)  :=  J  1  (5.4) 

The  arc  length  may  be  used  to  parametrize  points  on  the  center  line  of  the  cunent  configuration.  For 
the  convective  description,  this  is  not  necessary,  but  it  is  convenient  in  the  spatial  description.  We 
shall  tacitly  assume  that  the  rod  does  not  self-intersect;  that  is,  that  on  the  image  of  the  above 
mapping,  then-  is  a  well  defined  smooth  inverse  mapping  X  *-♦  S,  where  X  is  a  point  in  the 
image  of  <t>.  Tie  image  curve  is  parametrized  by  the  arc  length  and  so  we  will  regard  the  role  of  X 
as  being  played  by  S  in  what  follows. 


Motions  and  Velocity  Fields 

A  motion  is  a  curve  of  configurations  t  €  [0.T]  *-*  =  (<Pt.At)  e  C,  for  some  time 

interval  10, T].  Associated  with  a  motion,  one  has  the  material  velocity  field  V$(S,t) 
defined  by 

v*(S.t)  :=  |f  ♦(S.t)  :=  (<J>t(S.  t).  At(S.  t))  (5.4) 

Thus  at  any  time  t,  the  material  velocity  is  an  element  of  T^C,  the  tangent  space  to  Cat  the 


** i rw \  »v> 
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configuration  ♦  at  time  t.  The  spatial  velocity  field  V$(s,t)  is  defined  by 

V*(s,t):=  ($(S,  t).  A(S,  t)),  (5.5b) 

where,  as  above,  S  =  s(S,t)  is  the  arc  length  in  the  current  configuration.  Finally,  the 
convected  velocity  field  l>$(S,t)  is  defined  by  the  expression 

iVS.t)  :=  ([A(S.t)]T<(>(S.  t).  [A(S,t)]TA(S,  t))  (5.5c) 

Since  A(S.  t)e  S0(3).  we  can  write 

A(S.  t)*  A(S,  t)  W(S,  t)  =  W(s,  t)  A(S.  t)  (5.6) 

where  W(S,  t)  e  So(3)  and  W(s,  t)  6  So(3)  .  Accordingly,  the  convected  and  spatia 
velocity  fields  can  be  expressed  as 

IVS.t)  =  (  [A(S.t)]T  (f(S.  t).  *(S.  0) 

and  Vs.t).($(S.  t).  W($,  t))  (5.7) 

Remark.  One  can  easily  check  that  the  components  of  the  spatial  velocity  in  the  moving  frame 
(t,(S)}  are  equal  to  the  components  of  the  convected  velocity  in  the  inertial  frame  { E t } .  We  also 
remark  that  the  material  velocity  field  Vf  may  be  viewed  either  as  a  left  extension  of  the  convected 
velocity  field  or  as  a  right  extension  of  the  spatial  velocity  field  V*. 

Strain  Measures. 

As  in  the  three  dimensional  theory,  one  defines  the  (one-dimensional)  deformation  gradient 

as 


♦  ■(S.t)  :=  OP'(S.t).  A'(S.t)),  where  (  )'  :=  d/dS  (5.8) 

Similarly,  in  parallel  with  the  definition  of  velocity  fields,  starting  from  the  (Lagrangian) 
deformation  gradient  ♦  '(S.t)  one  defines  convected  and  spatial  strains  according  to 


rc 


* 


t 
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convected 

material 

spatial 

r(s.t) 

:=  [A(S,t)]T*<P‘  (S.t) 

<P‘(  S.t) 

2f(s.t)  :=  d<f>(S,t)/ds 

D(S.t) 

:=  [A(S.t)]T*A'(S.t) 

A‘(S.t) 

U(s,t)  :=  [dA(S,t)/ds]*AT(S.t) 

where  S  =  s(S,  t).  We  again  note  that  the  components  of  the  convected  strains  in  the  material 
frame  {E,J  coincide  with  the  components  of  the  spatial  strains  in  the  moving  frame  { t,(S) )  up  to 
a  factor  (because  of  the  arc  length  parametrization  that  may  be  used  in  the  spatial  description.)  The 
above  expressions  can  be  derived  from  the  three  dimensional  theory  by  a  duality  argument 
employing  the  formula  for  the  stress  power  given  in  the  remark  below. 


The  Equations  of  Motion  in  Spatial  Description 

Associated  with  the  motion  t  ►-*  4>t  e  Cone  assumes  the  existence  of  smooth  vector  fields 
n(s.t).  m(s,t).  and  a  scalar  field  p(s,t)  interpreted  respectively  as  the  contact  resultant  force, 
contact  resultant  couple,  and  density  per  unit  of  current  length.  These  fields  satisfy  the  following 
spatial  local  form  of  the  equations  of  motion: 


p.(jj-i)fLso 

PAv  =  ♦  n*  (5.10) 

•  ,  am  w 

p(jw  ♦  w  *  jw]  =  — y  ♦  o  *  n  ♦  m 

In  these  equations,  we  use  the  following  notation: 

d  8  8s 

5t  5s  "5t 

denotes  the  material  time  derivative,  and  7T ,  m  are  the  body  force  and  torque  per  unit  of  current 
lengths  =  S(S,  t).  The  role  of  the  Jacobian  is  played  by  J  -ds/dS;  J|(s.t)  is  the  time 
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t 

t 

/ 

i 

|  dependent  inertia  dyadic  of  the  cross  section  given ,  relative  to  the  moving  frame  ( t  j }  by 

J  *=  J*9  ♦  J33t3  ®t3  (5.11) 

j  where  J°^  and  J33  in  (5.1 1)  and  A  in  (5.10)  are  inertia  constants  associated  with  the  reference 

cross  section.  (Typically  :=  /q  dQ.  ot,0  =  1.2.  and  A  :=  Jq  dQ.  where 
£ 2  are  coordinates  in  the  reference  configuration  defined  by  (5.1).) 

Equations  (5.10)  are  the  local  statement  in  the  spatial  description  of  balance  of  mass,  balance  of 
linear  momentum  and  balance  of  angular  momentum,  respectively.  We  refer  to  Amman 
11972),[1976]  for  a  derivation  of  analogous  equations  from  the  three  dimensional  theory. 


Equations  of  Motion  in  the  Convective  Description 


In  the  convective  description,  the  balance  laws  are  expressed  directly  in  the  inertial  frame  { E  j } . 
To  this  end  one  defines  vector  fields  N(S,t),  M(S,t),  and  P„j(S)  as  pull-backs  of  their  spatia 
counterparts  by  means  of  the  relations 

N(S,t)  :=  [A(S,t)]T  -n(s.t). 

M(S.t)  :=  [A(S,t)]T  -mCs.t),  (5.12) 

P^S)  :=  J(S,  t)p(S.t) 

where  S  =  s(S,  t)  and  J(S,  t)  =  ds(S,  t)/dS  .  Omitting  explicit  indication  of  the 
arguments,  by  making  use  of  the  relation 

^  =  J-1  AIM'  ♦  Q  ■  M)  (5.13) 


and  the  analogous  expression  connecting  H  and  N,  a  straight  forward  calculation  then  yields  the 
following  statements  of  balance  of  mass,  linear  and  angular  momentum: 


dP^f 
dt  S 

Aprtfn>  ♦  w  *  v]  = 


0 


(5.14) 
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pJjW  ♦W*JW]=  ♦Q*M*r*N*M 

where  V  =  AT$  is  the  convective  velocity  ind  J  :=  ATjA  is  the  time  independent  inertia 
dyadic. 

Remark:  Within  the  context  of  the  three  dimensional  theory,  n,  and  m  are  defined  as  the 
resultant  force  and  the  resultant  torque  relative  to  the  line  of  centroids  of  the  distribution  of  stress 
acting  on  a  cross  section.  The  definition  of  2 f  and  6)  is  unique  in  the  sense  that  the  stress  power  is 
given  by 

Uuf ,f  dr  dS  *  ds  5  4>.uN'f  * M'61  ds 

(5.15a) 

Here,  P  is  the  first  Piola-Kirchhoff  stress  tensor,  F  is  the  deformation  gradient  of  the 
configuration  given  by 

2 

4>  =  <p(S,t)  ♦  I  £ata(S.t).  and  v=  '  -  W  *  (5.15b) 

<*=1 

is  the  co- rotated  rate  measuring  the  rate  of  change  relative  to  the  moving  frame. 

Appropriate  stress  measures  conjugate  to  the  Lagrangian  strains  (<P‘.  A')  may  also  be  obtained 
by  cither  (i)  left  extension  of  the  spatial  stress  measure  (n,m)  or  (ii)  right  extension  of  the 
convected  stress  measures  (N.M).  For  instance,  the  latter  extension  takes  the  form 

;=  Al?  and  :=  Afl  (5.16) 

Obviously,  such  an  extension  preserves  the  stress  power.  ■ 

Constitutive  Equations 

In  view  of  the  expression  (5.15)  for  the  stress  power,  one  characterizes  (isothermal) 
hypcrelastic  response  by  assuming  the  existence  of  a  free  energy  function  \ft(S,2f,(j),A)  such  that 


d^(S.y.CJ.A) 
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_  dplS.Z.O.A) 

"4  m  s  - eu - • 

By  postulating  that  these  equations  are  frame  indifferent;  i.e., covariant  under  the  left  action  of  the 
Euclidean  group  of  spatial  isometries,  a  further  reduction  is  obtained  as  follows.  Let  t  0(t)  e 
S0(3)  be  an  arbitrary  superposed  rigid  body  rotation,  and  let  $t4  :=  (  Q(t)<P.  Q(t)At)  €  C. 
Denoting  by  (*)*  objects  associated  with  4,  one  has  the  relations 

Z*  -  0Ct)2T,  n4  =  0(t)n 

6)4  =  Q(t)G>,  m4  =  0(t)m  (5.18) 

V  =  P 

Thus,  by  choosing  0  =  AT  ,  and  since  y/4(S,  02T,  Qcd,  QA)  =  S,  Z.  Li,  A),  it 

follows  that  the  only  possible  form  compatible  with  frame  indifference  and  locality  (see,  for 

example,  Marsden  and  Hughes  [1983])  is  given  by 


d'HS.r.Q) 

er 


e^(s.r.Q) 

8Q 


(5.19) 


Example.  An  example  of  a  constitutive  equation  consistent  with  the  above  invariance 
requirements  and  useful  in  computation  (see  for  example,  Simo  and  Vu-Quoc  [1986b])  is  furnished 
by  the  uncoupled  linear  systems  : 

N  =  C*(r  -  r°)  ;  M  =  Ch(Q  -  Q°)  (5.20) 

where  CN  and  are  symmetric  positive  definite  and  constant,  I*°  =  AqT<Po  wd  Q°  = 
AoTAq‘.  This  ensures  that  the  reference  configuration  $q(S)  ;r  (<P0(S),  Aq(S))  is  stress 
free.  Linear  constitutive  models  of  the  type  (5.20)  are  analogous  to  the  Saint  Venant  Kirchhoff 
model  of  three  dimensional  elasticity  and  are  typically  restricted  to  small  strains  due  to  the  behavior 
of  p.  We  refer  to  Ciarlet  [1986,  Chapter  4]  for  a  summary  of  appropriate  growth  conditions  of  p. 


Hamiltonian  in  the  convective  description 

One  defines  corvee  ted  linear  and  angular  momenta  according  to  the  following  expressions 


M  ^1*  ^  -  J  ^  , 
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M  ;=  PrefAV  and  T7  :=  p^fJW 


(5.21) 


Thus,  in  the  absence  of  body  forces  and  couples,  and  assuming  for  simplicity  pure  displacement 
boundary  conditions,  the  Hamiltonian  is  given  by 

H((r.  Q):  (n.  tt))  = 


]l(ca  p^-’SrtJ*  •  n-(  p^)-'n)DS 
0 

♦  J%(s,  r.  Q)ds 


(5.22) 


In  the  next  section  we  shall  consider  the  Hamiltonian  structure  of  this  rod  model  on  the  reduced 
space  T~  T*£/S0(3). 


"  ’  V  •  •  j-  _•  J 
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$6.  The  Hamiltonian  Structure  for  the  Geometrically 
Exact  Rod  Model  in  the  Convective  Picture 

In  this  section  we  shall  derive  the  Poisson  structure  which  makes  the  equations  for  a 
geometrically  exact  rod  model  Hamiltonian.  To  do  this  we  shall  proceed  in  a  way  that  is  similar  to 
that  in  three  dimensional  elasticity  and  the  free  rigid  body,  namely,  we  shall  reduce  the  canonical 
bracket  in  material  representation  by  means  of  spatial  isometries.  To  carry  this  out,  begin  by 
recalling  that  the  configuration  space  for  the  rod  mode!  is  given  by 

C=  {♦  |  ♦  =  (<P.A)  :  [O.L]  — *  R3  *  S0(3)  }  (6.1) 

The  tangent  space  to  C  at  a  configuration  ♦  is  given  by 

T tC=  (  V*  |  V*  =  <V*.  WA)  :  10.U  — • >  R3  *  TaS0(3)  }  (6.2) 

and  the  cotangent  space  is  given  by 

T iC=  {  M,  =  (M^p.  ltA)  :  CO.L!  —  R3  *  TaSO(3)  }  (6.3) 

Now  introduce  the  duality  pairing 

<Y>  :  T ,C*  T R 

defined  as  follows: 

<  V*.  M,>  =  | L  [  tyv,  •  J  i1a;*a]  is  («■«) 

Recall  that  the  reason  for  the  factor  1/2  is  that  according  to  (4.4)  and  (4.5), 

rtA:ftA  =  2  n  w 

Next,  we  turn  our  attention  to  the  formulation  of  the  canonical  Hamiltonian  structure  on  T  C 


(6  5) 
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The  Canonical  Bracket 


We  consider  functions  f.  g  :  C - *  R  so  that  the  canonical  Poisson  bracket  is  given 


Explicitly,  one  has  the  expression 


I'  .)  .(.M-S 


Sf  Sg 
6<J>  '  SMg, 


iL .  JL1  ds 

8<f> 


1  fLT  Sf  Sg  Sg  ST  1 

♦  -  J  I  YJT  :  "  u  :  ^ — J  dS  •  (6-7) 

2  Jo 1 5A  sftA  8A  sffA 

As  in  section  3,  we  remark  that  the  Poisson  bracket  (6.7)  is  written  as  if  the  variables 
and  (A,  ltA)  were  both  independent.  This  however,  requires  some  caution  since  lmC 
is  not  simply  a  produxt  space  Again,  we  refer  to  Abraham  and  Marsden  [1978]. 


Reduced  (Poisson)  Bracket 

In  the  present  context,  reduction  by  material  frame  indifference  amounts  to  considering 
functions  f  :  1*  C - ♦  R  which  are  S0(3)  invariant,  these  are  of  the  following  form. 

■  r((<p.  A);  (ty  rtA>) 

:  T((At»'.  ATA'):  (AtH#.  ATt1A)>  .  (68) 

where  (•)'  ■  d()/dS.  Introducing  the  notation 


T  :=  AT<p' 


fi  :=  ATA 
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n  :=  ATh0 

ll  := 


(6.9) 


which  is  consistent  with  the  left  representation  of  we  may  rewrite  equation  (6.8)  as 


f((<p.  A):  (M,p.  nA)>  =  f  «r.  C):  (rt.  ft)) 

=  T((r,  d)-.  (n.  n»  .  (6.io) 

One  now  obtains  a  reduced  bracket  in  terms  of  the  variables  ((1*.  Q);  {ft,  ft))  by 
transforming  the  canonical  bracket  with  the  aid  of  the  chain  rule.  The  key  result  is  contained  in  the 
following: 


Proposition  6.1.  We  have  the  following  formulae: 


8f 

8<p 


8f 

8A 


$T  er 

SM^  "  A  8n  ; 


dT 

“■sr 4 


n « 


dT 
dn  4 


l  _  dT  /  dT  \t 
s n  “  '5TT  *  (  iff) 


(6.11a) 


(6.11b) 

(6.11c) 


♦  Q  x 


(6.1  Id) 


Proof.  The  first  three  formulae  above  follow  directly  from  the  chain  rule.  To  prove  the  last  one. 
we  proceed  as  follows.  Making  use  of  the  chain  rule  we  have 


.  6A^  =JL{  —  *  6At<P‘  ♦ 

'  Jo  1  dr 


\  tf[  (5ATA-  ♦  AT8A  )] 


tTWVV mn  vi  ■v’  i. 
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♦  Ui '  tAT*‘  *  ?,rt  |fSATnA}dS 


(6.12) 


■  /0LATsA=[r®  4f«n®|±--ifi  |l]«.s 

•  V  ? SA:  IS-)  * A*  s  J  dS  •  <6 13) 


where  we  have  used  skew-symmetry  of  df/dfi,  integration  by  parts,  and  the  relation  A'  = 
A  Q  .  We  now  observe  that  since  SA  e  Ty^S0(3)  one  has  the  left  representation  5A  = 
ASS.  where  e  S0(3)  and  so  we  may  rewrite  (6.13)  as 


<  £  •  -> 


Jo  dr  dfi 

.in*  i(*L’] dS  . 
2  art  2  vefi'  J 


(6.14) 


However,  since  €  So(3),  we  have 


8§  :  A  s  S@  :  [ 


A  -  A 


for  any  A  e  GL(R3).  In  addition,  we  recall  the  identity 

[a  ®  b  -  b  ®  a]  =  •  [a  «  b]  A  ,  a,  b  e  R3  . 


(6.15) 


(6.16) 


along  with  the  relation 


[A6  -  6A3  =  IA  *  B]a  .  A.6e  so(3) 


Therefore,  since  d  I  7efi  =  IdT/dQf  ,  we  have 


(6.17) 
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<rk  SA>  *  -T  sS'Hr‘J!r‘n* 


dT  _  df 
Sn  '  °  *  ao 


1  „  d7  /  d7\  *1A 

4  2  n  *  i  dS  • 

2  dTT  dQ  (6.18) 


and  replacing  s$  by  AT6  A.  the  result  follows. 


Using  the  formulae  in  Proposition  6.1  and  employing  the  notation  { f .  g}  =  (f,  g)^  ♦ 
(f,  g)A  where 


(6.19a) 

(6.19b) 

we  have 


In  addition,  since  A  :  6/2  ■  A  ■  B, 


. .  ,  i  fL  /  df  dg  dg  df  \  .. 

{f- 9}a  :s  *  :artT)dS  • 

If.  g)»  °=  Jo  '  afi^'  dS  • 


9)a:= 


e7 
ST  4 


(6.21) 


f 


6 


I 


so  that 
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I*.  9)a  = 


/>•( 


dT  ^  dg  \  f  dT 

an  an 


>♦<&> 


„  /  d7  ag 

♦n'\  SS-*  AH 

♦  Q  •  ( 


eg 


d7  8  g  dg 


aQ  an  aQ 

4  r  (  dT  K  .ii  .  ii 

r  '  v  er  an  ar 

Thus,  we  have  proved  the  following 


dg 

an 

v  dQ 

an  ' 

an' 

*alf)] 

dS 

Theorem  6.2.  The  reduced  Poisson  bracket  on  IP- T*C/SO(3)  is  given  by 


(T,  g)  = 


(canonical) 


( interaction  ) 


flfD  ril.-tf-  il«  ill 

■  J„  Ur  ot  ar  an 1 


(Lie -Poisson  for  a 
semi-direct  product } 
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Note  that  the  first  integral  gives  the  canonical  term  in  the  variables  [(I*.  Q);  (M,  TT)];  the 
second  integral  gives  interaction  terms ;  and  the  last  rvwo  integrals  are  the  Ue-Poisson  terms  . 


Corollary  6.3.  Hamilton's  equations  f  -  {  f,  H  }  with  Hamiltonian  given  by  equation  (5.22) 
and  the  Poisson  bracket  given  by  (6.23)  are  equivalent  to  the  following  convened  equations  of 
motion 


where 


n  =N'*Q*N-w*n 

ri  sM'oQ*M*r*N-W*TT 

t  *  v*q«i;  .v/*r 
Q  S  W*  4  Q  X  w 


N 


er 


and 


(6.24) 


The  last  two  equations  in  (6.24)  can  be  checked  by  a  direct  calculation  from  the  kinematic 
relations  in  section  5.  The  first  two  equations  in  (6.24)  coincide  with  (5.14). 

In  summary,  we  have  found  the  reduced  space  !P-T*C/SO(3)  to  be  the  space  of  the 
convected  variables  ((I\  Q);  (FI,  TC))  with  the  Poisson  bracket  given  by  (6.23).  This  reduced 
bracket  has  been  obtained  from  the  canonical  bracket  in  material  representation  by  reduction.  As 
with  three  dimensional  elasticity,  if  the  motion  on  the  reduced  space  is  known,  then  the  original 
motion  in  the  material  description  can  be  reconstructed  by  solving  the  following  system 

_  AM  -  Al» 

at  "  AW  at  ’  Ay 


Remark.  The  Hamiltonian  formulation  can  be  of  assistance  in  the  establishment  of  conservation 
laws,  in  addition  to  being  useful  for  stability  and  bifurcation  studies.  For  example,  suppose  that  the 
rod  has  an  isotropic  cross  section,  so  that  the  inertia  dyadic  J  has  two  equal  eigenvalues,  just  as  in 
the  case  of  the  Lagrange  top.  Then  there  is  a  material  symmetry  group  acting:  it  is  the  group  S? 
acting  on  the  right .  Namely,  let  the  symmetry  group  S1  consist  of  rotations  about  the  symmetry 
axis,  say  U.  The  group  action  is  then  given  by  the  following  action  of  R  e  S*:  and 

unchanged,  and  (A,  flA)  -  (AR,  flAR).  Here  the  momentum  map  corresponding  to  this 
action  gives  the  conserved  quantity 


J7.  Geometrically  Exact  Plate  Model 


In  this  section  we  summame  the  basu  equations  go.em.ng  the  goemetncally  exact  plate 
model.  First,  we  introduce  some  necessary  notation 


Notation 

In  the  geometric  description  of  the  configuranons  of  a  plate,  the  unit  sphere 

S2  :=  {t  6  R3  |  |t|  =  1}  .  (7.1a) 

plays  a  central  role.  Given  any  te  S2,  the  tangent  space  at  t  is 

TtS2  =  (Vt  €  R3  |  Vt  ‘  t  =  0}  .  (7.1b) 

Thus,  TtS2  can  be  identified  with  R2.  In  addition  to  S2  one  introduces  the  following  subset 

of  S0( 3)  which  enables  one  to  establish  a  link  between  the  plate  and  rod  formulations. 

Let  E  €  R3  be  a  fixed  but  otherwise  arbitrary  vector.  Define  SE  to  be  the  set  of  rotations 
A  6  S0(3)  whose  rotation  axis  e  R3  is  normal  to  E;  i.e., 

SE  :=  (A  €  S0(3)  |  y  e  R3  satisfies  A^  =  ^  and  >jt  •  E  =  0}.  (7.2) 

The  set  SE  is  closely  related  to  S ,  as  the  following  (well-known)  result  shows. 

i 

Proposition  7.1.  Given  any  t  e  S  there  is  one  and  only  one  element  At  e  SE  such  that 

t  =  AtE  .  (7.3a) 

In  fact,  S2  and  SE  are  diffeomorphic. 

Proof.  Given  t  €  S  define  the  orthonormal  set  {Ejl  by 
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E2  5=  Es  *  E,  . 

It  follows  that  Ej  =  t  -  (E  *  t)E.  Clearly,  {Ej}  is  uniquely  defined  for  t  *  E.  Let  ©  := 
COS^tE  •  t)  and  define  €  S0(3)  by 

At  :=  E)  $  E)  «  cos  ©(Ej  ®  E2  ♦  Ej  ®  Ej) 

♦  sin  ©(E2  ®  Ej  -  Ej  ®  E2)  • 

Hence,  At  €  S£ .  Uniqueness  follows  from  the  construction.  ■ 

The  geometric  interpretation  of  this  proposition  should  be  clear.  It  constructs  At  by  rotating 
E  to  t  in  the  plane  they  span,  through  the  angle  0  . 

The  tangent  space  to  SE  at  the  identity  is  given  by 

T,Se  ;=  {§  e  So{3)  |  0  *  E  =  0}  .  (7.3b) 

Finally,  the  tangent  space  TASE  at  A  €  SE  is  obtained  by  using  either  left  or  right  tranlsadons 
in  SO(3).  For  instance,  for  left  invariant  vector  fields  we  have 

TASE  =  {£a  :=  A@  j  forany  TjSe,  and  Ac  SE}.  (7.3c) 

An  analogous  characterization  holds  for  right  invariant  fields. 

Kinematic  description  of  the  plate 

We  consider  the  kinematic  description  of  a  plate  with  thickness  h  >  0.  Essentially,  no 
conceptual  modification  is  required  for  the  more  general  case  of  a  shell.  Let  {Ej(X°),  1=1, 2. 3} 
be  a  fixed  orthonormal  basis  for  R3,  with  E3  ■  E.  Any  point  X  e  R3  may  be  expressed  as 

X  =  X°  ♦  U,  where  X°  =  X-E*.  £,  €  R  .  (7.4a) 

From  the  point  of  view  of  the  three  dimensional  theory,  the  reference  configuration  of  the  plate  is 
these;  SC  R3  given  by 
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*  '*  { x  =  X°  •  tE  I  X°  «  D  and  £  c  [-£.£]  }  ,  (7.4b) 

where  Q  C  R2  is  a  given  region  in  the  plane.  One  refers  to  {X°  ♦  £E  |  X°  €  dQ  and  £, 
c  [-  h/2  .  h/2]  }  as  the  edge  of  the  plate,  and  to  Q  as  its  mid-plane  .  The  basic  kinematic 
assumption  is  that  any  admissible  configuration  $  s  3  — *  R3  of  the  body  is  characterized  as 

X  =  ?(X)  :=  <J>(X°)  ♦  £t(X°)  .  .  (7.5a) 

where  <P  :  Q  — »  R3  maps  the  mid -plane  onto  <P(Q)  €  R3,  and  t  is  a  unit  vector  attached 
to  the  point  <P(X°)  €  R3,  not  necessarily  normal  to  <P(Q),  which  is  referred  to  as  director  or 
fiber  direction.  The  fact  that  t  need  not  be  normal  takes  into  account  possible  shear 
deformations. 

Note  that  the  mapping  X°  e  0  — •  t(X°)  6  S2  assigns  to  points  X°,  vectors  t  in  the 
unit  sphere  (at  <p(X0)).  By  Proposition  7.1,  the  unit  vector  t  at  <P(X°)  can  be  obtained  from 
E  uniquely  through  a  rotation  A(X°)  c  Sg  with  rotation  axis  normal  to  E. 


Configuration  space 

The  kinematic  assumption  (7.3)  embodies  two  essentia]  ingredients:  (i)  Points  X°  €  Q  in 
the  midplane  arc  mapped  onto  points  X  €  R3  through  the  mapping  <P  :  Q  — *  R3;  and(ii)unu 
vectors  E  €  S2  attached  to  points  X°  e  Q  are  mapped  into  unit  vectors  t (X °)  e  S2  to 
points  <P(X°)  e  R3  by  rotations  A  :  Q  — »  Sg  vith  rotation  axis  normal  to  E.  According 
to  this  view,  two  abstract  characterizations  of  the  set  C  of  possible  configurations  of  the  plate  arc 
possible: 

(a)  Director  point  of  view: 


C  :=  {  ♦  ■  (<p.  t)  |  <p  ;  D  — *  R3  and  t  *  Q  —»  S2  }  .  (7.6a) 

The  tangent  space  T^-  C  at  a  configuration  ♦  €  C  is  defined  as 

c  :=  {  6?  .  (8<P,  6t)  |  5(p  :  Q  —  R3  and  5t  :  Q  —  TtS2}  . 
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(b)  Constrained  frame  point  of  view.  Equivalently,  as  a  result  of  Proposition  7.1,  the 
configuration  space  C  may  be  defined  as 

C  :=  {$  ■  (<p,  A)  |  <P  :  Q  — *  R3  and  A  !  Q  — »  SE}  .  (7.7a) 

According  to  this  view,  the  tangent  space  is 

T $c  :=  (s'*  a>  (Sep.  5A)  I  S<P  :  Q  — *  R3  and  SA  :  Q  — ♦  TAS£}.  (7.7b) 

Points  of  view  (a)  and  (b),  although  equivalent,  lead  to  different  forms  of  the  equations  of 
motion.  Starting  with  the  classical  paper  of  Ericksen  and  Truesdell  [1958],  a  substantial  body  of 
contemporary  work,  see,  e.g.,  Amman  [1976,1978],  Green,  Naghdi,  and  Wainwright  [1965],  and 
Naghdi  [1972,  1980]  has  been  typically  concerned  with  the  director  point  of  view  (a).  HeTe,  on  the 
other  hand,  we  take  a  different  approach  and  adopt  the  constrained  frame  point  of  view  (b)  as  a 
starting  point.  Our  motivation  for  this  lies  in  the  strucure  of  the  equations  of  motion  and  the 
Poisson  bracket  in  approach  (b)  and  which  are  essentially  identical  to  that  of  the  geometrically  exact 
rod  model  considered  in  Sections  5  and  6.  The  form  of  the  Poisson  bracket  corresponding  to 
approach  (b),  in  the  convected  description,  is  obtained  by  a  reduction  process  that  amounts  to 
enforcing  explicitly  the  additional  (symmetry)  condition  that  the  rotation  fields  must  lie  in  the  subset 
SE  C  S0(3). 


Motion  and  velocity  fields 


Amotion  is  a  curve  of  configurations  t  €  [0,  T]  — ►  =  (<Pf,  Af)  e  C  for  some 

time  interval  [0,  T]  C  R4.  Associated  with  the  motion,  one  has  the  material  velocity  field 
V$(X°,  t)  defined  as  usual  as 

V*(X°.  t)  :=  Jr  *(X°.  t)  ■  («j>t(X°),  At(X0))  .  (7.8a) 


Thus  the  material  velocity  is  a  mappting  t  »-*  V 


€  T $  C  where  T is  the  tangent  space 


at  configuration  ♦.  In  addition  we  define  the  convected  velocity  field  l*$(X0,  t)  e 
*(Q)  as  in  section  5,  by  the  expression 
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V*(X°.  t)  :=  ([A(X°.  t)]T<i>t(X°,  t).  [A(X°.  t)]TA(X°,  t))  .  (7.8b) 

Here,  X(Q)  denotes  the  space  of  smooth  fields  :  Q  — »  R3  *  TjSe.  Since  A(X°,  t) 

€  SE.  we  have  A(X°.  t)  =  A(X°,  t)W(X°,  t)  where  ft(X°,  t)  £  TtSE.  Thus,  the 
convected  velocity  field  can  be  expressed  as 

1VX°.  t)  =  (V,  W)(X°.  t)  (7.8c) 

where 

V(X°.  t)  :=  [A(X°.  t)]TV(X°.  t)  .  V(X°,  t)  :=  t)  . 

Note  that  the  material  velocity  field  V$  is  the  left  extension  of  the  convected  veclocity  field 
Vq.  Similarly,  one  could  define  a  spatial  field  whose  second  factor  w(X°,  t)  is  defined  as 

W(X°,  t)  :=  A(X°,  t>[A(X°,  t)]T  .  (7.8d) 

Then,  W  :  Q  *  [0,  T]  — »  TjSt  is  interpreted  as  a  spatial  angular  velocity  ,  where 
TjSt  :=  {0  €  so(3)  |  t  '  0  =  0}  is  the  tangent  space  at  t  =  AE. 

Strain  measures 

As  in  the  geometrically  exact  rod  model  outlined  in  Section  5,  starting  from  the  LagTangian 
deformation  gradient  ♦  ^(x0,  t),  (erf  =  1.  2) .  one  defines  convected  strains  according  to  the 
expressions 


Convected 

TJX*.\)  :=  At<P  *(X°,  t) 


Material  Spatial 

<P.a(X».  t)  V" 


Oa(X°.  t)  =  [VAJX0.  t)]v 


A  *(X°.  t)  =  [A.^T  ]V 


We  note  that  the  above  expressions  can  also  be  derived  from  the  three  dimensional  theory  by  a 
duality  argument  employing  equivalence  of  the  stress  power.  This  is  discussed  briefly  in  a  remark 
below  (see  equation  (7.15). 
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Stress  resultants  and  stress  couples ; 

equations  of  motion  in  spatial  description 

Associated  with  the  motion  t  *-*  €  C,  one  assumes  the  existence  of  smooth  vector 

fields  n^X0,  t),  m^X0.  t),  (o<  =  1.  2),  and  Pref(X0.  t)  interpreted  as  internal  resultant 
force,  internal  resultant  torque,  and  density  per  unit  of  area.  These  fields  satisfy  the  following 
spatial  local  form  of  the  equations  of  motion 


* n  =  p^v  , 


"V*  * 


<p 


(ot  *  no<  4  m 


=  PrttKW 


(7.10) 


where  n".  m  are  the  applied  body  force  and  torque,  and  Prefh,  prefK  are  intertia  coefficients. 
Typically,  for  constant  thickness  plates.  Pref K  =  prelh3/12.  Note  that  W’t  :  0.  that  is  W  e 
TjSt  for  each  X°  e  O  and  t  €  10.  T].  Equations  (7.10)  are  the  local  statement  in  the  spatial 
description  of  balance  of  linear  momentum  and  balance  of  angular  momentum,  respectively.  We 
refer  to  Libai  and  Simmonds  [1983]  for  a  derivation  of  thpse  equations  from  the  three  dimensional 
theory. 


Remark.  From  the  point  of  view  of  the  three-dimensional  theory,  the  right-hand-side  of  (7.10; 
agrees  with  the  standard  definition  of  linear  and  angular  momentum  per  unit  of  reference  surface 
relative  to  the  mid-plane.  To  see  this,  consider  configurations  $:$--*  R3  of  the  form  <|>  ;= 
<p  ♦  £,AE,  where  E,  e  [-  h/2  ,  h/2].  The  angular  momentum  TC  relative  to  the  mid-plane 

is  then  given  by 


7t 


r  h/2  , 

(#  -  <P)  *  P„f<P  ^ 

-h/2 


r  h/2 
-h/2 


pre£AE  «  ($  ♦  AE) 


=  -^~a[e  ■  atAe] 
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Pi.!*3 

12 


ME  ■  WE) 


12 


AIE  *  (W 


E)3 


P^3 

12 


AW 


p^’ 


(7.11) 


where  we  have  used  the  fact  that  We  T  |Sg  so  that  W  *  E  =  0 .  ■ 


Convective  equations  of  motion 

An  alternative  form  of  the  balance  laws,  which  is  particularly  convenient  in  computational 
applications,  is  obtained  from  the  convective  description.  As  far  as  we  are  aware,  this  form  of  the 
equations  has  not  been  stated  explicitly  in  the  literature.  As  in  the  rod  model  discussed  previously, 
the  basic  idea  is  to  pull-back  the  spatial  balance  laws  (7.10)  with  the  orthogonal  transformation  A 
:  Q  — *  Sg.  Accordingly,  one  defines  vector  fields  N(X°,  t),  M(X°,  t)  by  the  expressions 


N(X°.  t)  :=  ATn(X°.  t)  ,  M(X°,  t)  :=  ATm(X°,  t)  .  (7.12) 

By  making  use  of  the  relation 

mot,J  =  ASM*,  ♦  *  MJ  ,  (7.13) 

and  the  analogous  expression  relating  and  N^,  a  straightforward  calculation  yields  the 
following  statements  of  linear  and  angular  momentum: 

p„,h[l>  •  W  »  l>]  =  Na  oI  .  Q*  «  N„  .  N  .  (7.14a) 

P^KW=  •  ra  «  Na  •  M  .  (7.14b) 


Remark  .  Within  the  context  of  the  three  dimensional  theory,  and  represent  resultant 
forces  and  resultant  torques  relative  to  the  mid -surface  of  the  distribution  of  stress  acting  on 
sections  S^.  As  in  the  rod  model  discussed  above,  and  are  uniquely 

determined  in  the  sense  that  the  stress  power  is  given  by 
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JQ„j.  »  «,  p  :  *  d0  =  JQ  ["«  “«]  OQ 

f  ,  <7  15) 

=  JQ  [N„  '  f«  *  Ma  •  Oa]  dQ  . 

where  P  is  the  first  Piola-Kirchoff  stress  tensor,  F  the  deformation  gradient  of  the  configuration 
<P  :=  <P(X°,  t)  ♦  £»t ( X t),  and  (  )v  =  {  )'  -  W  *  (')  is  a  co-rotattd  rale.  ■ 


Further  reduction  of  the  convective  equations  of  motion 

In  contrast  with  the  rod  model  considered  in  Section  5,  the  convective  equations  of  motion 
(7. 14)  aur  amenable  to  further  reduction.  The  reason  for  this  additional  reduction  --  which  can  be 
earned  out  either  in  the  spatial  or  in  the  convective  descriptions  —  is  that  only  Sg  C  SO(3) 
enters  in  the  configuration  space  of  the  plate,  not  the  entire  S0(  3)  as  in  the  rod  model. 

Proceeding  in  the  convected  description,  to  carry  out  the  reduction  of  equations  (7.14)  we 
first  note  that  since  Q^,  €  Tj$£,  it  follows  that  Q^*  E  s  E  ■  0  and  so  * 

is  parallel  to  E.  We  exploit  this  fact  by  introducing  the  decomposition 

K  n0°*o0e  .  r*  =:  r’-zj..  (7.16) 

A  straightforward  calculation  then  yields 

ra -  n_  =  rjm  h>.  (r  » .  E)oa -  (nJ>. e)e0  , 

(7.17) 

Q*  «  N*  =  Q*  «  Nj  ♦  (Qot  *  E)°*  • 

Making  use  of  (7.16)  and  (7.17),  balance  of  angular  momentum  (7.14b)  now  reduces  to 

Qo<  *  Mo<  4  *£  *  N£  =  0  ’  (7.18a) 

•  (rM°-  E)oa  -  (n,>  e)e„  .  M  »  =  p„^W  .  (7.18b) 

The  structure  of  equations  (7.18)  suggests  the  introduction  of  the  following  notation: 
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SO 


w° 

:=  O^nE 

:=  W  «  E 

or 

o* 

w 

II  II 

m  m 

*  . 

« w°  , 

«2 

:=  M^-E 

Mot 

ii 

m 

«M®  . 

In  addition,  we  denote  by  P£  :=  I  -  E  $  E  the  onhogal  projection  parallel  to  E.  Since  W  * 
V  z  (  E  *  W°  )  *  V,  the  convected  acceleration  may  be  expressed  as 

AT  #  r  [l>°  ♦  ^W°]  ♦  [  V-  W°  V°]e  ,  (7.20a) 


where 


V°  :=  P^V  ,  and  V  :=  E  ■  V  .  (7.20b) 

We  then  have  the  following  five  equations  of  balance  of  momentum  in  the  convected  representation: 

<«,  *  QX,  •  N  «  =  Prrfh[l>°  .  VW°  ]  . 

■  K  *  0  =  Pld*[>-  W°  t>°]  .  (7.21) 

Hi.-  -  r«°o„  •  Nj>=,  •  M  0  =  PrdKW »  . 

Remark.  Interpretation  of  W  °and  Q^.  The  equations  of  motion  (7.21)  along  with 
conditions  (7.18a)  are  in  fact  the  convected  form  of  the  equations  of  motion  corresponding  to  a 
director  description  of  the  plate,  as  in  (7. 6a,b).  The  vector  field  W°  is  the  convected  director 
velocity  as  the  following  identity  shows: 

ATt  =  AtAE  =  WE  ■  W  *  E  =:  W°  . 


An  entirely  analogous  interpretation  holds  for  Q^.  ■ 
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Constitutive  equations 


In  view  of  expression  (7.15)  for  the  stress  poweT,  one  characterizes  (isothermal)  hyperelastic 
response  by  assuming  the  existence  of  a  free  energy  function  of  the  form 


wx0.  rSt.  =*,  tf*). 

As  in  the  rod  model  considered  in  Section  5,  this  form  of  the  free  energy  in  the  convective 
description  arises  naturally  by  postulating  invariance  under  left  action  of  the  Euclidean  group;  i.e., 
material  frame  indifference.  We  then  have  the  equations 


.  av(x  rg°.  g,.  oa) 


(7.22a) 


dv(x°.  r°.  s„.  q,°) 


(7.22b) 


dv(x  °.  r°.  5*.  Q^) 


(7.22c) 


We  note  that  equation  (7.18a),  arising  from  balance  of  angular  momentum  about  E  can  be 
expressed  as 

Qj. »  «  n,°  =  o. 

From  this  condition  one  obtains  the  following  restriction  on  the  constitutive  equations  (7.22): 

Q0„  av(x».rjsM.oj)  r„  dv(xMt,.so,,0°)  =  g 


(7.23) 


Condition  (7.23)  is  analogous  to  the  symmetry  of  the  Cauchy  and  second  Piola-Kirchoff  stress 
tensor  in  three  dimensional  continuum  mechanics,  which  in  tum  results  from  balance  of  angular 


momentum. 
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Hamiltonian  in  the  convective  picture 


One  defines  the  convected  linear  and  angular  momenta  by  the  expressions 

n  =  n°  +  m e  :=  prefMi>°  ♦  ve] 

TT  :=  preJKW 

n°  ::  n  K  E  -  PrefKW°  . 


(7.24) 


Thus,  in  absence  of  body  forces  and  assuming  for  simplicity  homogeneous  Dirichlet  boundary 
conditions,  the  Hamiltonian  is  given  by 


H  ■■  5  JQ  {<  P^KlM2  •*■]•<  PjO-’In*!*}  on 


(7  25) 


Convective  Equations  of  Motion  in  Director  Notation 

For  completeness  we  record  below  the  field  equations  governing  the  proposed  plate  model  in 
director  notation.  As  in  Section  7.2,  let  {Ej}  be  the  standard  basis  in  R3,with  Ej  s  E,  {X^} 
covering  Q  C  R2,  and  Xj  ■  £  e  [-h/2,  h/2j.  Define 


m°  :«  M  °  ®  E*  ,  :=  N°  ®  E*  ,  ©  :=  0^  , 

ID  0  :=  Q°  ®  E*  ,  ®°  :=  r£  ®  E*  ,  X  :=  5 JE*  , 


(7.26) 


The  momentum  equations  (7.21)  may  then  be  recast  in  the  following  form; 

vnJl0  ♦  ID0©  ♦  H°  =  prefMl>°  ♦  VW°]  , 

VIV fi  -  1D°  :  Ji°  ♦  o'  =  p  fhli'-  W°  •  V°]  ,  (7.27) 
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<DIVM°  -  ©°©  ♦  ft°X  ♦  M°  =  prefK W°  . 


Xhis  director  form  of  the  equations  is  particularly  convenient  for  the  treatment  of  constitutive 
equations,  as  the  following  example  illustrates. 

Example  A  particular  class  of  constitutive  equations  satisfying  the  invariance  requirement  (7.23) 
may  be  constructed  by  an  argument  that  mimics  that  of  the  three  dimensional  theory.  First,  we 
assume  the  following  more  restrictive  condition  that  implies  (7.23): 


0,5 «  =  0  .  and  r°  «  N.J  =  0  .  (7.28) 

(a)  Constitutive  equations  for  N5, .  Assume  that  ^  depends  only  on 
Define  the  following  stress  tensors: 

:=  R°  ,  :s  R°(e°)T  .  (7.29) 

Then,  condition  (7. 28)^  is  equivalent  to  the  symmetry  of  and  ^1°;  that  is, 

ra°  »  N.J  *  0  ~  *°  =  (8-V  »  ®°  =  (®°)T  (7.30) 

-r  0  "ITT'  ^ 

To  complete  the  development,  we  define  strain  tensors  conjugate  to  and  sL  by 

®°  ==  i  [(®°)T  ®°  - 1]  .  [i  -  (®°)'T®°  -’] 

(7.31) 

A  direct  calculation  shows  that  these  tensors  satisfy  the  stress  poweT  relation 

0  0  0  ■  0  0  ■  0 
N  o<  •  r*  =  $  :  £  z  Z  :J)  .  (7.32) 


Accordingly,  we  define  constitutive  equations  of  the  form 


evex0, «  ) 
e®° 


z{ 


e^{x°.  I ?°> 
e*>° 


(7.33) 
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The  net  result  is  an  uncoupled  set  of  constitutive  equations  that  satisfy  at  the  outset  restriction 
(7.23).  In  particular,  one  may  consider  a  linear  relation  between  W  and  ©  . 

(b)  Constitutive  equations  for  M^.  It  is  possible  to  construct  an  uncoupled  equation  for 

of  the  form  (7.22b)  satisfying  at  the  outset  the  restriction  (7.28)j  .In  particular,  this  relation 


can  be  linear. 


§8.  The  Hamiltonian  Structure  for  the 
Geometrically  Exact  Plate  Model 


In  this  section  we  develop  the  Hamiltonian  structure  for  the  plate  model  in  the  convective 
representation  with  governing  equations  summarised  in  the  previous  section  .  We  show  that  these 
equations  are  Hamiltonian  relative  to  a  non -canonical  Poisson  structure  in  the  space  2*  of  convected 
variables 


{a-0*.  e#i  q»„).  (  n°.  m,  n°)  ). 

The  First  Reduced  Bracket 

Our  derivation  of  the  corresponding  Poisson  bracket  follows  the  same  reduction  scheme 
employed  both  in  three  dimcndional  elasticity  and  in  our  treatment  of  geometrically  exact  rods,  and 
can  be  outlined  as  follows.  We  start  with  the  canonical  Hamiltonian  structure  on  the  cotangent 
bundle  of  the  configuration  space,  T *C,  and  the  corresponding  canonical  bracket.  From  this 
bracket  one  obtains  by  left  reduction  by  SO(3);  i.e.,  by  material  frame  indiference,  a  Poisson 
bracket  for  the  plate  model,  which  is  the  counterpart  of  the  Poisson  bracket  derived  for  the  rod 
model.  An  important  difference,  however,  is  that  one  can  further  reduce  this  bracket  by  enforcing 
the  condition  that  only  a  pan  of  SO(3),  namely  Sg  C  SO(3),  enters  in  the  configuration  space. 
This  additional  reduction  is  the  result  of  an  extra  symmetry  of  the  Hamiltonian,  which  is  now 
invariant  with  respect  to  rotations  about  the  axis  t.  Physically  ,  this  symmetry  corresponds  to  the 
fact  that  no  inertia  and  no  stiffness  is  associated  with  rotations  about  the  director  t . 

The  bracket  in  terms  of  the  variables  {(I*,*,  Q^),  (7*1,  TT) }  is  derived  from  the  canonical 
bracket  in  material  variables  in  the  same  way  as  for  the  rod  model.  This  bracket  is  as  follows: 


f  \(  d7>) 

l  *9 

(  dg 

\  dT 

|T.  9) 

,c an 

'oc '  an 

(canonical) 

sg  t 

■35^V 

SR")  dQ 

(8.1a) 

(interaction) 

■y»- 1  £ 

„il. 

an 

dg 

(8.1b) 
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•J„(°.  t 


(Lie-Poisson  for 
a  semi-direct 

product) 


dg 

dg 

M  dT  1 

W  " 

dQot 

dn  J 

dg 

dg 

«  dT  1 

en 

dr  * 

dn  J 

d7  d  g  dg  dT 


rdf  dg 
♦  n  ■  Lam  *  5n  ' 


an  an  an  an 


-L*-lUr 


87  _  dg 


an  “  5tT 


]  dQ 


]}  10 


(8.1c) 


The  Further  Reduced  Bracket 

The  further  reduction  is  acomplished  by  using  the  change  of  variables  in  (7.19)  given  by 


r*  =  r0**  =  je  n  =  n°  ♦  me 

q*  z  e  *  q0*  n  =  e  *  TT° 

Since  £•.&  e  TjSg  and  ft  e  TjSg,  we  have  the  constraints 
TT  '  E  =  0  and  D^’  E  =  0, 

and  therefore 


Q*  E  «  Q°ot 


(8.2a) 

(8.2b) 


•  ll 
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From  (8.2a)  and  the  chain  rule,  we  get 


ef  er  et_ 
s  ar°a  ‘  an. 


(8.6) 


where  f  is  a  function  of  {(I^,  Q^),  (M,  TT)  ),  as  in  the  bracket  above,  and  f  is  a  function  of 
the  variables  {(r®^.  He<t  (  M®,  M,  T7®)  )  which  are  related  to  the  preceding  variables 

through  equations  (8.2).  Similarly  one  has 


8g  _  djf 

err0  "  an 


(8.7) 


ef  _  er  E  ar 
dn  an0  dM 


(8.8) 


Substitution  of  (8.4)  -  (8.8)  into  (8.1)  gives  the  following  bracket: 


Theorem  8.1.  The  reduced  bracket  on  the  space  T  of  the  variables 

{cr0^.  q°j.  ( n°.  m.  n«) } 


is  given  by 

(f.g)  - 


U°V 

L  en° 

dg 

dS* 

dg 

en° 

er  i 

■ShJ 

♦  a  V  1 

r  df 

dg 

dg 

er 

*  rV  1 

■  df 

■  ^ 

eg 

en° 

dg 

er 

en° 

(interaction) 
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(canonical) 


♦  =  .  f  df 

.  dg 

dg 

.  df 

01  L  an0 

to 

an0 

er0* 

.  J  {n"  -  [  4t 

JD 

eg 

dg 

ar 

an0 

dlW0 

en° 

(Lie  Poisson  for  a 

semidirect 

♦  *•[  V 

dg 

*g  . 

a,l 

product) 

an0 

an0 

an0 

dn°  J 

]} 


.  f  f(  *!_ )  .^s_.(is_)  ,*L 

lr\  iv£q0^/,o<  dTT°  V  ■*tt0 


'Q 


dQ°  •*  dTT° 


.  ( i!_)  .  *£. .  (  a9  )  .  1L 

~  N-n  "  en° 


drv“  dn°  dr°„  •“ 


.  (  .iL)  .  *£. .  (ii.)  .  ii  }  d! 

'  '  nl  S  s%  S  W  J 


as*  •*  a*r  as*  ■* 


As  for  the  rod  model  discussed  in  Sections  5  and  6,  from  the  expression  (7.25)  for  the 
Famihonian  and  the  bracket  (8.9),  we  have  the  following: 

Corollary  8.2.  Hamilton's  equations  in  the  form  f  =  {  f,  H  }  ,  where  f  is  an  arbitrary 
function  on  the  phase  space  given  by  theorem  8.1,  with  the  bracket  given  by  equation  (8.9)  and 
Hamiltonian  given  by  equation  (725)  are  equivalent  to  the  following  convective  equations  of 
motion: 


*f-Q o,  *-  q^h\  ♦  w°n° 

Tt.0  -  u0  4n0  q  .  anjO 
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v«'  w»-ra-QVv°  (8.10) 

.  wo=a .  Q^v 

where  we  have  set  N0^  =  d^/dl*0^,  0^=  d^p/dE^  and  M0^  =  d^/dD0^. 


The  first  three  of  these  equations  of  motion  have  been  derived  already  (see  7.21).  The  second 
group  of  three  equations  can  be  directly  checked  using  the  kinematic  relations  given  in  section 
7.The  corollary  then  follows  by  the  general  principles  of  reduction,  but  may  also  be  verified  by  a 
direct  computation. 


Conclusions 

In  this  paper  we  have  presented  a  systematic  development  of  the  Hamiltonian  structure  for 
geometrically  exact  nonlinear  elasticity,  including  solids,  rods,  and  plates.  We  have  emphasised  the 
convective  representation  since  it  is  numerically  convenient  and  it  is  useful  for  the  coupling  with 
rigid  body  dynamics  to  obtain  models  for  rigid  bodies  with  flexible  attachments.  In  particular,  our 
formulation  of  the  geometrically  exact  plate  model  in  the  convected  representation  constitutes  the 
natural  counterpart  of  the  classical  Kirchhoff-Love  and  Reissner-Antman  models  for  rods.  The 
derivation  of  the  Poisson  structure  follows  the  same  lines  as  earlier  works  of  Marsden,  Ratiu, 
Weinstein  and  their  coworkers,  namely  it  is  obtained  by  reduction  from  canonical  brackets  in  the 
material  (or  Lagrangian)  representation  by  reduction;  ie,  by  the  elimination  of  rotational  symmetry 
(material  frame  indifference),  and  by  introducing  the  Cauchy-Green  tensor  as  a  dynamical  variable. 
This  approach  is  consistent  with  the  covariance  investigations  of  nonlinear  elasticity  by  Hughes, 
Marsden,  and  Simo.  In  a  future  work  we  shall  use  the  Hamiltonian  structures  for  geometrically 
exact  rods  and  plates  to  study  the  nonlinear  stability  of  rigid  bodies  with  flexible  attachments, 
following  the  ideas  of  Krishnaprasad  and  Marsden. 

It  is  clear  from  the  literature  and  related  work  that  these  methods  are  much  mote  significant 
and  widely  applicable  than  the  results  given  here  may  directly  indicate.  For  example,  the  attention  to 
the  proper  geometry  and  the  nonlinear  context  that  is  typified  by  the  present  investigation,  is  of 
benefit  for  numerical  work,  as  has  been  demonstrated  by  Simo  and  VuQuoc.  Also,  other  models 
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can  be  investigated;  for  example,  it  is  clear  from  the  literature  (Iwinslri  and  Turski.  Marsden  and 
Weinstein)  that  one  can  also  include  electromagnetic  effects  into  the  same  formalism. 


%  '  '  ■  A.W  Wl A  M 


•"*  ir»  y» rrjm y- y-«J/»VF_,  *>  <JVJ1W  V V  v  v.v.  r; 
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